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Abstract 

.^ ' The average number of minijets and the corresponding transverse energy pro- 

cn ■ duced in heavy ion colhsions are evaluated by including explicitly semi-hard parton 

^jj ! rescatterings in the dynamics of the interaction. At the LHC semi-hard rescatterings 

(^ I have a sizable effect on global characteristics of the typical inelastic event. An in- 

^^ ■ teresting feature is that the dependence on the cutoff which separates soft and hard 

^-^ parton interactions becomes less critical after taking rescatterings into account. 
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1 Introduction 

Semi-hard physics is one of the most important issues in the interaction of heavy ions at the 
LHC Given the total energy involved and the large number of participants, the component 
of the inelastic interaction which can be described within a perturbative approach is in fact 
rather substantial in heavy ion collisions at the LHC |l], 0. The result is the production 
of a large number of minijets in the typical inelastic event. The description of the semi- 
hard component of the interaction, adopted by the majority of the papers on the subject, 
follows the approach of ref . : the semi-hard component of the inelastic cross section, in 
a collision of two nuclei with atomic mass numbers A and B, is written as: 

^ j j.p(''.,T.Bm\^-.,T,m , (1.1) 



n=l' 



where 






lxx's>Apf^ rr, 

TABiP) = ld%rA{h)rB{h-l3). (1.2) 



In Eq. (|1.2| ) r^(&) and Tsip — [3) are the thickness functions of the two interacting nuclei, 
normalized to A and B, respectively, and depending on the transverse coordinates of the 
interacting partons h and b — (3, where (3 is the impact parameter of the nuclear collision. 
aj is the single scattering cross section to produce jets, expressed as a convolution of 
the nuclear parton distributions G^(a;), Gb{x') and of the partonic cross section a-^-l" 
(integrated on the momentum transfer with the cutoff po and without the multiplicity 
factor of the produced jets). The indices / and /' label the different kinds of interacting 
partons and the momentum fractions x, x' are defined with respect to the single nucleon 
momentum. The dependence on the scale factor is implicit in all quantities. 

The semi-hard cross section, as expressed in Eg . ( |1 . 1]) , represents the contribution 
to the total inelastic cross section of all events with at least one semi-hard partonic 
interaction, so that one may write ainei = o'n + o'soft, the two components being separated 
by means of Pq. Eq.( p..l| ) may be derived in an eikonal approach to nucleus- nucleus 
interactions, by writing the eikonal phase as the sum of a soft and of a hard component 
|1|, Q. The physical picture corresponding to Eq.( |l.l| ) is that of a distribution of multiple 
independent parton collisions localized in different points in transverse space and with 
the average number depending on the nuclear impact parameter. The average number of 
parton interactions, at fixed P, is (iV(/3)) = ajTAsiP) so that, if aj is at the lowest order 
in as, the average multiplicity of minijets in a given nuclear collision (that is at (3 fixed) is 
2{N{/3)). The integrated inclusive cross section for producing minijets is therefore given 
by 2 J d^l3{N{P)). To obtain the inclusive cross section one needs, in fact, to count the 
jets produced with their multiplicity, so that, in the case of multiple parton interactions, 
the inclusive cross section is obtained by taking the average of the distribution in the 
number of parton collisions. The result is that the integrated inclusive cross section is 



given by the single scattering expression, cjj, multiplied by the average multiplicity of jets 
produced in a single partonic interaction, which shows that the description of the process, 
as given by Eg. ( p. .11 ), is consistent with the AGK cancellation [^ (all unitarity corrections 
cancel in the inclusive). The cancellation property obviously holds for all averages, that 
are therefore equal to the result obtained by means of the single scattering expression, so 
that the transverse energy produced by minijets is given by 

{E,m = 2TABmj PtP^d\. (1.3) 

Jpt>PO " Pt 

The semi-hard cross section an is a smooth function of po for small values of the cut-off: 
the limiting value of an is in fact the geometrical limit ti^Ra + RbY, Ra and Rb being 
the two nuclear radii. On the contrary {N{j3)) and {Et{j3)) are singular at small Pq and 
their behaviour may be roughly estimated on dimensional grounds to be {N{l3)) ~ 1/pq 
and {Et{(3)) ~ 1/po- The singular behavior at low pq can be used to set the limits of 
validity of the picture. Indeed, for the picture of the interaction to be valid one should 
take a relatively large value of the cut-off po! in this way the whole semi-hard interaction 
takes place in a relatively dilute system and the overall number of interactions will be 
relatively small. To deal with a regime where the number of parton interactions and 
the density of the interacting partons are large, the main modification adopted by the 
majority of papers is to include shadowing corrections in the nuclear parton distributions 
1^, 0. In this way one obtains a substantial reduction of the number of projectile and 
target partons at low x and the picture can be extended to sizably lower values of the 
cutoff po- Even so, when po is further reduced, one reaches the condition of a highly dense 
interacting system where the whole picture ceases to be valid, which sets the lower limit 
for a sensible choice of the cutoff po [|, ^ . 

The overall resulting features are therefore that po is different when varying the atomic 
mass number of the interacting nuclei and their energy, and that the distribution in the 
number of hard collisions at a fixed value of the impact parameter /3 is a Poissonian in 
the whole semi-hard regime so that all average quantities are computed, as above, with 
the single scattering expression. 

The clean physical interpretation of the approach, which incorporates the geometrical 
features of the nuclear process, unitarity, the factorization of the hard component of 
the interaction and the AGK cancellation rule, justifies the great success of the picture. 
Still there are a few delicate points which deserve further investigation and where the 
description of the process might be improved. 

A very general approach to nucleus-nucleus interactions is trough the Reggeon Dia- 
gram Technique, where at high energies the interaction is described by the exchange of 
many Pomerons, including both independent exchanges between different nucleons and 
multi-Pomeron interactions, which represent the collision of a given projectile nucleon 
with several different target nucleons in a given interaction process. Each partonic colli- 
sion corresponds to a fluctuation with a large transverse momentum inside the structure 
of an exchanged Pomeron. While the simplest case is that with a single partonic loop 
with large p^, whose discontinuity corresponds to a 2 — > 2 partonic collision, one might 
imagine more complicated fluctuations, with several connected parton lines, all with large 
Pt, in the structure of a multi-Pomeron interaction. The simplest case of this kind is 
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Figure 1: a) A triple pomeron interaction between 3 nucleons. b) Parton ladders with the topology of 
Fig.nk). The thick lines are fluctuations with high pt and the dashed line is a possible cut. c) Partonic 
view of the cut of Fig.0b): P is a projectile parton interacting with two target partons T^ and the final 
state is given by the three large-pt jets J^. 



illustrated in Fig.|I|, where the Pomerons are represented by parton ladders. Discontinu- 
ities of such fluctuations originate conflgurations, where the compensation of transverse 
momenta in the flnal state involves several large pt partons, while the whole large pt con- 
flguration is generated by partons belonging to different chains, representing Pomerons 
attached to different nucleon lines. In the picture of the semi-hard interaction just re- 
called those conflgurations contribute to the shadowing corrections of the nuclear parton 
distributions. Their discontinuities, on the contrary, are not included in the semi-hard 
interaction dynamics. On the other hand, when approaching the black disk limit, the 
initial and flnal partonic states become locally isotropic in transverse space, so that both 
initial and flnal state conflgurations need to be built up with lots of partons. A proper 
discussion of semi-hard dynamics in the black disk limit requires therefore taking into 
account fluctuations with many large pt parton lines, which are to be included not only 
in the virtual corrections but also in the production process. Hence one needs to con- 
sider partonic interactions where several partons, with low virtuality and pt and sizable 
longitudinal components (so that, in the nucleus-nucleus cm. frame, each one may be 
ascribed to a deflnite nucleus) interact producing large pt jets. The simplest possibility 
of this kind (one projectile against two targets, t/s -^ and basically electro-dynamical 
interaction between point like objects) was discussed in Ref . [p!0| , where the forward am- 
plitude of the process and all the cuts were derived. The flnal result of the analysis was 
that, in the t/s —>■ limit, the different cuts of the amplitude are all proportional one 
to another, the proportionality factors being the AGK weights. It was moreover possible 
to express the three body interaction as a product of two body interaction probabilities. 
Therefore, one may argue that, in the black disk limit and in the nucleus-nucleus cm. 
frame, a projectile parton interacts with all the various nearby partons of the target in 
the different directions in transverse space, and that the corresponding cross section may 
be approximated by the product of two-body interaction cross sections, so that one might 
call such a process re-interaction or re-scattering. 

An attempt to introduce such more elaborate semi-hard dynamics, including explicitly 
semi-hard parton rescatterings in the interaction, was done in ref. |jll| and [1^ . Both the 



average number of minijets and the average transverse energy are modified by semi-hard 
rescatterings, and an interesting feature is that both quantities develop a less singular 
dependence on the cutoff, in such a way that the choice of po becomes less critical when 
semi-hard parton rescatterings are taken into account. The average number of minijets 
and the transverse energy produced in heavy ion collisions have been recently discussed 
in several papers, with the purpose to determine the initial conditions for the further 
evolution and termalization of the system (see e.g. 0, rB|). We think that it might 



be interesting to have an indication of the effects of rescatterings on these quantities at 
the LHC, and we'll discuss this topic in the present paper. In the next section, after 
including rescatterings in the picture of the interaction, we'll recall the expression of the 
average number of minijets and derive the corresponding average transverse energy. Then, 
we'll give some quantitative indication on the effect and comment the qualitative features 
induced by the more structured interaction dynamics. 



2 Multiple Parton Scatterings and Average Number 
of Minijets 

The introduction of semi-hard parton rescatterings can be obtained in a picture of the 
interaction where the soft part is factorized in a Poissonian multi-parton distribution and 
the hard part is expressed in terms of perturbative parton-parton collisions. The nuclear 
multi-parton distributions are unknown quantities and the reason of the choice of the 
Poisson distribution is that it corresponds to the case where the information on the initial 
state is minimal, since the whole distribution is expressed in terms of its average value 
only. One has moreover the possibility of introducing systematically further informations 
on the nuclear partonic structure in terms of correlations among partons [Q . The hard 
part is written in terms of two-body collisions by introducing the probability of having at 
least one interaction between the two configuration, in a way analogous to the expression 



of the inelastic nucleus-nucleus cross section ||T5[. The process is therefore represented as 



the sum of all possible interactions between all configurations with a definite number of 
partons of the two nuclei. In this way the semi-hard component of the cross section can 
be written as 



^H = fd'pf:^T^{x,M-P)...rAx^,b^-P)e-^d'^d%r^{x,b-(3) 



n=l 

oo 



E \ r^K, h\) . . . r^(x;, &;) e" / dx'^b'Tsix', v) 



=1 

n I 



where 



1 — ]^ ]^(1 — aij) dxid%i . . . dxnd'^bndx[d'^b[ . . . dx'id^b'i , (2.1) 

j=i i=i 



TA{x,b)=TA{b)G{x) 



To keep the notation as simple as possible, the indices labeling the different kinds of 
partons have been suppressed and the dependence on the cutoff po is implicit. The two 
Poissonian distributions, with average numbers TA_{x,b — P) and TB{x',b'), represent the 
multi-parton distributions of the two interacting nuclei. The probability to have at least 
one partonic interaction, given a configuration with n and / partons in the two nuclei, is 
represented by the square parenthesis in Eq.( p.l|) and, according with the discussion in 
the previous paragraph, is constructed by means of the probability atj of interaction of 
a given pair of partons i and j. The interaction probability is a function of the cutoff 
Po, so that only the interactions with momentum exchange larger than po are taken into 
account in (|2.1|) . The majority of interactions takes place with a momentum exchange 
close to the cutoff value, we therefore evaluate the parton distributions at the scale of the 
cutoff. Since the distance over which the hard interactions are localized is much smaller 
than the typical nuclear radius, the probability of interaction can be approximated by 
d-{xiXj, hi - bj) = a{xiXj)5^'^\hi - bj). 

The most important features of Eq. ( p.l|) are that all possible interactions between the two 
partonic configurations are included, and that probability conservation is explicitly taken 
into account by the term in square parentheses. The Eikonal cross-section in Eq.(|l.l|) is 



obtained when one neglects all rescatterings in Eq.( |2.1| ), in such a way that each parton 



is allowed to interact at most once [|I6|. If, on the contrary, one keeps rescatterings into 
account one cannot write a closed expression for an- However, it is possible to obtain 
simple expressions from Eq.( p.l| ) for many relevant quantities. 

If one works out the average number of parton collisions {N{f3)) one obtains, as in the 
traditional approach, the single scattering expression (iV(/3)) = ajTAsiP) [Q- So the 
overall average number of parton collisions satisfies the AGK cancellation and is not 
affected by any unitarity correction; however it is not a simple quantity to observe. A 
more directly observable quantity, or at least one which can be more directly related to 
observable quantities, is the average number of produced minijets. An important effect of 
including rescatterings is that the number of produced minijets is no more proportional 
to the number of collisions, because now each projectile parton is allowed to interact more 
than once with the target. As a consequence, while unitarity corrections do not change 
the average number of collisions, they affect the average number of minijets produced in 
the nuclear collision. 



The correction term can be derived in a straightforward way from Eq. (|2.1|) ||11|| , but 
one may use also a more direct argument. One may in fact obtain the usual semi-hard 
cross section, Eq.( |l.l|) , by starting from the single scattering cross section to produce 
large-pi jets. In the perturbative QCD-parton model the cross section is written as: 



aj= (f(3d^bdxdx'TA{x,b-P)a{xx')TB{x',b)= d^p{N{p)) . (2.2) 

Jxx's>ipo J 

The expression needs however to be unitarized also when the cutoff po has a rather large 
value, since aj is proportional to the large factor A x B. Eq. ( |1 . 1| ) is in fact the result of 
the s-channel unitarization of aj and may be obtained by noticing that in Eq.( p.2D (A^(/3)) 
is dimensionless and may be understood as the average number of parton interactions at 
a given impact parameter j3. The s-channel unitarized cross section an is the result of 
replacing this average number with the interaction probability, which, if the distribution in 
the number of interactions is a Poissonian, is just l — exp[—{N[(3))). Hence the unitarized 



cross section an represents the contribution to the total cross section of all events with 
at least one couple of partons interacting with a transverse momentum exchange above 
Po, as it is clear from the second line of Eq.(|l.lD. On the other hand aj, that includes 
also the multiplicity of the interactions, represents rather the integrated inclusive cross 
section (apart from the factor representing the average multiplicity of jets produced in a 
single collision). 

When the cutoff is moved towards low values and rescatterings need to be taken into 
account, the average number of jets produced is no more proportional to the average 
number of collisions. In this case one may proceed by applying to {N{f3)) an argument 
analogous to that previously used to unitarize aj. By looking at Eq.(|2.2|) one can identify 



(jiB{x,h)) = / dx Tb{x ,h) a{xx') (2.3) 

J xx' syApQ 

as the average number of collisions of each interacting A-parton at fixed x and h. Then 
one can write the average number of produced minijets at fixed impact parameter as 



2{N{I3)) = / d%dxTA{x, b - 13) {neix, b)) + / d%dx'TB{x', b) (^^(x', 6 - /?)) , 

(2.4) 
that represents the average number of incoming partons from the nucleus A multiplied 
by the average number of collisions against the partons of nucleus B, plus the analogous 
term with A and B interchanged. Then, if one replaces in Eq. ( |2.4| ) the average number 
of scatterings of each parton with its interaction probability, viz. 1 — exp{{nB{x, b))) and 
1 — exp((n^(x, b — P))), one obtains the average number of "wounded partons" Wa{x, b) 
and Wb{x', b — (3) of the two nuclei. These ones are the partons of the two nuclei that had 
at least one hard interaction. The expression of the average number of wounded partons 
of nucleus A (with transverse coordinate b and momentum fraction x, in an event with 
nuclear impact parameter (3) is therefore 

WA{x,b,f3) = TA{x,b-l3)[l- e-<^^(^'^))] . (2.5) 

Every wounded parton obtained in this way produces a minijet in the final state and the 
transverse energy produced by semi-hard interactions is the transverse energy carried by 
the wounded partons. As a consequence both the average number of minijets and their 
average transverse energy are quantities affected by the presence of rescatterings and the 
corresponding correction term is more and more important when the average number of 
scatterings {hb) becomes larger and larger, namely at low values of the cutoff po and (or) 
for large atomic mass numbers. 

The overall number of produced minijets, i.e. the sum of the wounded partons of nucleus 
A with those of nucleus 5, obviously coincides with the usual result 2{N{(3)), when the 
average number of rescatterings is small. When the number of rescatterings is large the 
two quantities are however different and, while the average number of collisions (iV(/3)) 
may be divergent in the saturation limit, the average number of wounded partons is on 



the contrary well defined. In fact one obtains that the square parenthesis in Eq. (|2.5| ) has 
1 as a limiting value and, in this limit, the average number of wounded partons is just 
the sum of the average number of partons of the two interacting nuclei. 
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Figure 2: Average number (jib) of semi- hard scatterings suffered by a parton from the A- nucleus while 
it is interacting with the B-nucleus. We consider a Pb-Pb collision at a cm. energy ^/s — 6 TeV/A, with 
a scale Q = po/2 and a fc-factor k = 2.5; Woods-Saxon thickness functions are used. Top: {ub) is plotted 
as a function of the projectile parton fractional momentum x in the case of a central parton-nucleus 
collision. The three curves are the result for a cut-off po — 2 GeV (solid line), 3 GeV (long dashes), 4 
GeV (short dashes). The shaded area corresponds to the central pseudo-rapidity region \rj\ < .9, where 
T] = \og{x y^s) / pq) . The dotted line marks the transition value of 1 scattering per parton. Bottom: (ub) 
is plotted as a function of the impact parameter for a parton emerging with pseudo-rapidity 77 = 0. 



To have a quantitative feeling for the regime of interest at LHC energies we plot (n^), 
the average number of scatterings per parton from the A-nucleus, as a function of the 
projectile parton fractional momentum x and impact parameter b. As one can see, at a 
cut-off po = 2 GeV a parton in the central region experiences on average 1 to 3 scatterings 
over the whole target transverse area except the very peripheral regions. As discussed 
above, at higher cut-off values the number of scatterings decreases. As a result, semi-hard 
rescatterings should be negligible down to po — 3 GeV where their effects begin to show 
up, and already at po = 2 GeV they should have a large effect and should not be neglected. 
These conclusions will be confirmed by the numerical computations discussed in section 

i 

In summary, the effect of rescatterings on the average number of the produced mini- 
jets is to reduce the number obtained by means of the single scattering expression, not 
differently, qualitatively, from the result of including shadowing corrections in the nuclear 
parton distribution. On the other hand the overall distribution in the number of minijets 
produced is modified. In the traditional approach there is a strong correlation in the dis- 
tribution of the number of minijets, since each minijet has a recoiling companion; when 
the average number of rescatterings increases this correlation gets weaker and weaker, 
so that in the high density limit no correlation is left and the distribution tends to a 



Poissonian [O. A further important difference is that, since Wa and Wb are well defined 
in this limit, after including rescatterings the average number of minijets becomes much 
less dependent on the choice of the cutoff at low pq. 
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3 Average Transverse Energy 



As we have discussed in the last section, the average number of minijets is not the only 
quantity modified by this more elaborate interaction dynamics, which in fact has a non- 



trivial effect also on the transverse energy carried by the minijets [|12[. The wounded 
partons of nucleus A at a given x and b and in a nuclear interaction with impact parameter 
/3, Eqn.f pISD are obtained by multiplying the average number of partons of A (with given 
X and b) by the corresponding interaction probability, which is a Poisson probability 
distribution in the number of scatterings, with average {nB{x,b)) = J dxT b{x' ,b)a{xx'): 



WA{x,b,(3) = TA{x,b-(3)J2 



{nBix,b)) {nj,{xM 



u=l 



V\ 



(3.1) 



We may therefore obtain the differential distribution in the transverse momentum pt of the 
wounded A-partons by introducing a constraint in the transverse momentum intergrals 
that give the total cross sections in the above expression: 

dWA 



d^Pt 



{x,b,P) 



TA{x,b-(3)f:- fTs{x[,b)...TB{x'^,b)e-Idx'Ts{x',b)a{xx') . 



da 



da 



(5(2) (ki 



^u — Pt) d^ki . . . d^ky dx[ . . . dx'j^ . 



(3.2) 



d'^ki d^ky 

The limits of integration on x'^ and x' are respectively xx'^s > Akf and xx's > Ap^, and 
all the distribution functions are evaluated for simplicity at a fixed scale. To obtain the 
corresponding average transverse energy {pt{x,b, P))a one has to integrate Eq . (|3^) with 
an additional factor pt. A convenient way to proceed is to introduce the Fourier transform 
of the parton-parton scattering cross section 

"^"^ a(0) = a (3.3) 



aiu 



d ke 



jku 



d^k 



and to express the 5-function in Eq. (|3.2|) in an integral form. Since (t{u) depends only on 
the modulus of u one obtains 



/•OO /"OO 

{pt{x,b,(3))A = TA{x,b-(3) dptpl duuJoiptu) ■ 

Jo Jo 

^-Jdx'TB{x',b)a{0) 



(3.4) 



OO -| 

^ — / dxTsix', b)a{u) 



u=l 



ui 



After multiplying the integrand by the exponential factor exp(— Api) one can exchange 
the two integrals and first do the integral on pt. The result of this integration is 



f 

Jo 



-^pt 



MPtu)Ptdp 



2X' 



u 



1 d 



(A2+m2)5/2 y^^ 



(A2+m2)5/2 



(3.5) 



which allows one to perform the integral on u by parts. After differentiating the series in 
V one obtains the factor 



a'iu) 



/■OO rjfj 

-2n k^J,(ku)—dk, 
Jo d'^k 



(3.6) 



which is proportional to u because of the argument of the Ji Bessel function. The inte- 
grand in u is therefore regular for u = also in the limit A — ;> 0. Then one takes this limit 
and sums the series in z/. The final expression has a closed analytical form: 

{p,(x, b,P)) = - Ta{x, b-P) r du dx'~^^VB{x\ 6) e / d^"^B{x\ h) [d{u) - a(0)] ^37) 



u 

The average transverse energy in an event with a given impact parameter (3 is the result 
of integrating Eg . (13 .71 ) on h and x and of summing the two contributions of the wounded 
partons of the nuclei A and B. Notice that the expression in Eq.( p.7| ) is much less depen- 
dent on the choice of the cutoff po than the usual average energy evaluated with the single 
scattering expression of the perturbative QCD parton model: the Rutherford singularity 
of the parton-parton cross section is in fact smoothed in Eq. ( p.6|) by the Bessel function 
Ji{ku), in such a way that the dependence on the cutoff po is only logarithmic. The same 
logarithmic dependence on the cutoff is present in the argument of the exponential 

da 



a (u) - a (0) = 27V f[jo{ku)-l\-^kdk (3. 



as a consequence of the behavior of [Jo{ku) — 1] for k —* 0. 

In summary, when the semi-hard cross section is expressed by Eq.( p.l|) one may obtain, 
without further approximations, a closed analytic form both for the average number of 
minijets and for the corresponding average transverse energy and in both cases the singular 
dependence for small values of the cutoff is smoothened by rescatterings. 



4 Discussion 

Many papers have been recently devoted to the production of minijets in heavy ion col- 
lisions [0, 1^ H, |13[. A rather general feature is that, because of the singular behavior 
of the elementary parton interaction at low momentum transfer, many relevant quanti- 
ties depend rather strongly (typically like an inverse power) on the value of the cutoff 
which distinguishes soft and hard parton interactions. The feature is unpleasant since 
although one might find physical arguments to determine a meaningful value of the cutoff 



P, 1^, |T^, it is rather difficult to fix it in a very precise way. We have therefore tried, 
in the present paper, to face this issue by studying the effect of a more elaborate inter- 
action dynamics on the average number of minijets produced in a nuclear collision and 
on the corresponding average transverse energy. While in the traditional picture of the 
semi-hard processes each parton is allowed to interact with large momentum transfer only 
once, we have included semi-hard parton rescatterings in the dynamics of the interaction. 
Semi-hard rescatterings, that are negligible when the threshold between hard and soft 
processes is high, become more and more important when the threshold is lowered and 
the target approaches the black disk limit. Naively one would expect that the inclusion of 
rescatterings in the picture of the interaction might worsen the divergent behavior at low 
transverse momenta; on the contrary a more careful analysis, that takes probability con- 
servation consistently into account, shows that the result is just the opposite. Following 



111 and [0 we have in fact represented the semi-hard nuclear cross section with Eq. ( pTT] 
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where all possible multi-parton collisions, including rescatterings, are taken into account 
and the conservation of probability is explicitly implemented. The average number of 
minijets and the corresponding transverse energy, at fixed x, b and impact parameter /3, 
are then expressed in a closed analytic form by (|2.5| ) and (|3.7| ), whose behavior with the 
cutoff is much less singular in comparison with the analogous averages obtained without 
taking rescatterings into account. The reason of this smoother behavior is that rescatter- 
ings introduce (through Ta and Tb) a new dimensional quantity in ( P75| ) and ( |3.7| ), the 
nuclear radius, which gives the dimensionality to the two average quantities at small po. 
When rescatterings are neglected the dimensionality at small po is provided by the cutoff 
itself, and the result is that the two quantities behave as an inverse power of the cutoff 
for Po -^ 0. 

Apart from the approximation of writing all connected multi-parton processes as prod- 
ucts of 2 ^ 2 partonic collisions, our approach states on the assumption of neglecting 
production processes at the partonic level (namely 2 — > 3 etc. parton processes) and of 
using forward kinematics in the nucleon-nucleon cm. frame. To have a feeling on the 
validity of such approximations at LHC energy, we have evaluated the average energy of 
a partonic interaction in the parton-parton center of mass frame, and the average value 
of momentum fraction a; of a projectile parton: 



{Ec.m.)crj = / dxdxWxx'sy^^G^^{x)a^-^ {xx')Gf^{x') 

ff 

{x)aj = / dxdx'xy^^G/^{x)a^^ {xx')Gj^{x') (4.1) 



When the whole rapidity range is considered typical values are (Ecm.) ~ 25GeV, (x) ~ 
3 X 10^^ (corresponding to a momentum of ~ 45GeV, if the nucleon-nucleon cm. energy 
is 6TeV) with po = 2GeV. The average (x) becomes substantially smaller when averaging 
in a narrow window in the central rapidity region. The relatively low value of (Ecm.), 
as compared with the cutoff, indicates that the inclusion of 2 —>■ 3, or of higher order 
partonic processes, should not spoil the whole approach, that could therefore represent 
a reasonable lowest order approximation. The relatively large value of the momentum 
boost to go from the nucleon-nucleon cm. frame to the partonic cm. frame shows, on 
the other hand, that forward kinematics is reasonable in the former frame of reference. 

The effect of semi-hard parton rescatterings on the average number of minijets and 
on the average transverse energy produced in a central Pb-Pb collision at LHC energies 
is summarized in Fig.^. We plotted the average number of minijets and their transverse 
energy in the case of a very central rapidity window, \y\ < .9, corresponding to the 
ALICE detector, and in a larger rapidity interval, \y\ < 4, that will be covered by the 
CMS detector. Figure |^a) shows, in the two cases, the dependence of the average number 
of produced minijets on the choice of the cutoff po- The dashed curves are the results 
obtained by the single scattering expression, Eq. (|2.2|) , while the continuous curves are 
the result of the inclusion of semi-hard parton rescatterings, Eq.( p.5|) plus the analogous 
term for the 5-partons. These expressions have been computed with the GRV98LO 



distribution functions [|T8| with no shadowing corrections included, and by representing 
the elementary partonic interaction at the lowest order in QCD; to account for higher 
order corrections the result of the elementary interaction has been multiplied by a factor 
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Figure 3: Average number N of minijets, a) b), and average transverse energy Et, c) d), in a Pb-Pb 
central collision in two different rapidity windows, \y\ < 4 and \y\ < .9, as a function of the cutoff, a) and 
c), and of the fc-factor, b) and d). The dashed curves are computed without including rescatterings, the 
continuous curves after including rescatterings. 



k = 2.5. Both the fc-factor and the scale Q = Po/2 where chosen in order to reproduce 
the value of the pp mini-jet cross-section at ^/s = 900 GeV |jT9[. The sensitivity to the 
value of the fc-factor is shown in Fig. ^) , where the curves have the same meaning as in 
Fig. ^), and the cutoff has been fixed to the value po = 2 GeV. Analogous curves for 
the average transverse energy carried by the produced minijets are shown in Fig. ^) and 
^). The average transverse energy without rescatterings has been computed by using 
Eq.( |1.3|) (dashed curves) and with rescatterings (continuous curves) by using Eq.( p.7|) , 
after integrating on b and on x (inside the corresponding rapidity windows) and adding 
the analogous contribution of the 5-partons. 
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a) 



b) 



Figure 4: Production of minijets in the transverse momentum plane, a) The projectile parton (P) is 
allowed to interact only once with a target parton (T), so that the two produced minijets are back-to-back, 
b) The projectile can scatter against more than one target and minijets are produced with a star-shape. 



5 Conclusions 



In this paper we have discussed the production of minijets in heavy ion colhsions at 
the LHC, by including exphcitly semi-hard parton rescatterings in the dynamics of the 
interaction. The regime of interest is, in fact, the regime where the nuclear target reaches 
the black disk limit for a projectile parton interacting with momentum exchange above the 
cutoff po (which defines the lower limit in pt for minijet production). In the mechanism of 
production of minijets usually considered (a projectile parton interacts with only a single 
target parton) the elementary interaction generates a state with a preferred direction in 
transverse space (the direction of the two minijets) which is not consistent with local 
symmetry in the transverse directions implied by the black disk limit. The star-like shape 
in transverse space of the state produced by the multiple interactions of a projectile with 
different target partons (scattered projectile -|- recoils) recovers, on the other hand, the 
symmetry property of the interaction in the black disk limit, see Fig.§. 

The basic element in our estimate is to recognize that the usual expression that gives 
the average number, {N{i3)), of parton collisions at a fixed value of the nuclear impact 
paramete, Eq. (|2.2| ), is obtained by convoluting the average number of partons of the pro- 
jectile with the average number, {nsix, b)), of interactions of each projectile parton with 
the target nucleus, Eq. (|2.3| ). Notice that the evaluation of {nB{x,b)) and its dependence 
on the momentum fraction is, in this way, determined in a unique way by the parton distri- 
butions. Our results are therefore parameter free. Given the expression of the semi-hard 
cross section, Eq. (|2.1| ), the average number of produced minijets, W{x, b, /3), Eq.( p.5|) , and 
the corresponding average transverse energy {pt{x,b, j3)), (Eq.^.TD, are computed without 
further approximations, so that our result for A^ and Et, plotted in the figures, are exact 
consequences of the nuclear cross section ( |2.1| ). An approximation that is done when 
writing the nuclear cross section ( |2.1| ), is to evaluate the parton distributions at the scale 
of the cutoff, so that in Eq. ( |3.2|) the only dependence on the transverse momenta is in the 
elementary partonic cross sections, which basically corresponds to neglecting the logarith- 
mic dependence of the distributions in comparison with the inverse power dependence of 
the cross section. However, our formalism can be extended to take into account a general 
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scale dependence. 

The main features of the numerical evaluation are that semi-hard rescatterings have 
a sizable effect on the average number of minijets and on the transverse energy produced 
in heavy ion collisions at the LHC, so that they affect also global characteristics of the 
typical inelastic event. The induced correction increases with the value of the /c-factor, 
which represents higher orders in the elementary parton collision, and with the size of 
the rapidity window, since rescatterings are more frequent for partons with a larger mo- 
mentum fraction. By looking at the dependence on the cutoff, both the average number 
of minijets and the corresponding average transverse energy are more regular at low po, 
showing a tendency to saturate below po — 3 GeV and making in this way the choice of 
the cutoff less critical. 
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